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Abstract 

The blackbody radiation in a box L 3 with periodic boundary conditions in thermal equilibrium 
at a temperature T is affected by finite-size effects. These bring about modifications of the thermo- 
dynamic functions which can be expressed in a closed form in terms of the dimensionless parameter 
L T. For instance, when LT~4 - corresponding to the value where the most reliable SU(N) gauge 
lattice simulations have been performed above the deconfining temperature T c - the deviation of 
the free energy density from its thermodynamic limit is about 5%. This may account for almost 
half of the pressure deficit observed in lattice simulations at T ~ 4T C . 

PACS numbers: ll.10.Wx; 11.15.Ha; 12.38.Mh 
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In a very hot quark gluon plasma, when the temperature is much larger than any other 
relevant mass scale, asymptotic freedom leads to expect that the effective coupling to be 
used in thermodynamical calculations should be small. However, even in the case where 
the coupling is very small, strict perturbation theory cannot be used, the reason being that 
infrared divergences occur in high order calculations and various resummations are needed to 
get meaningful results. Lot of effort has been devoted to the perturbative calculations of the 

. The values obtained by adding successively high order contributions 
oscillate too much, and strongly depend on the renormalization scale. Thus such a plasma 
cannot be described simply as a gas of weakly interacting quarks and gluons. 

This very conclusion was also reached by lattice studies both for pure SU (3) gauge case 



{(J and for different kinds of fermions [7|, [8|, [9[ . Such calculations revealed a slow approach to 
the ideal gas limit of the thermodynamic functions. In particular it resulted a large deficit in 
the pressure and entropy as compared to the Stefan-Boltzmann law for free gluon gas, which 
remained at the level of more than 10% even at temperatures as high as T ~ 4T C . Similar 
results have also been found for £77(4) and ££7(8) gauge theories in a more restricted range 
of temperatures |XQf ] . 

These simulations were made on lattices of size N^xN t with periodic boundary conditions. 
Much effort has been dedicated to study and control the ultraviolet (UV) cut-off effects which 
are in general 0((aT) 2 ). In the standard Wilson formulation, temporal extent N t = 1/aT > 
8 is needed in order to get reliable extrapolations of the thermodynamic functions to the 
continuum limit. 

In this paper I wish to focus on another facet of lattice simulations, i.e. the infrared 
(IR) finite-size effects. In fact, for a thorough comparison of the numerical data of the 
hot quark-gluon plasma with the Stefan-Boltzmann (SB) law, one should consider a free 
gluon gas enclosed in a box with the same size and the same boundary conditions of the 
corresponding numerical experiment. It is clear that finite size effects are expected to be 
particularly relevant in a free boson gas: the lack of an intrinsic length scale leads to a 
maximal sensitivity to the geometrical shape N s /N t = LT of the system. The purpose of 
this paper is to evaluate these infrared effects. 

This might appear to be an academic exercise in view of the fact that in T — > oo limit the 
internal energy density e of the resulting ideal lattice gluon gas has been explicitly evaluated 
Uj through numerical integration both in the thermodynamic limit and for fixed N s /N t . 
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The finite-size effects turn out to be of the order of 1% for N s /N t = 4 [11] (see Fig.l). 
Moreover, in the continuum limit this system is scale invariant, therefore the trace of the 
energy-momentum tensor is vanishing 

T^ = e - 3p = (T -> oo, a -> 0, iV s /^ = const) , (1) 

hence the deviations of the pressure from the SB value are even smaller. 

There is however a missing point in the above reasoning. The quantity which is calculated 
in lattice simulations is not the pressure p, but the free energy density /, the reason [ljj] 
being that the evaluation of p would involve the derivative of the bare coupling with respect 
to the volume which is known only perturbatively, while / can be evaluated in a sounder way 
and in the thermodynamic limit one has p = —f. In a finite volume, however, this relation 



is violated. A numerical study on a SU 
suggested the empirical rule L s /L t > 4 



'2) gauge system at intermediate temperatures 
| to get rid of the finite-size corrections. In 



fact this bound is not enough at high temperature: in the free gluon gas limit, where the 
canonical partition function Z can be evaluated exactly even on finite volume, we shall prove 
that 

?r2 (LT) 3 — log VuT + 0(e~ 27TLT ) . (2) 



N 2 — 1 45 

As expected, this quantity is not purely extensive, owing to the finiteness of the volume 
V = L 3 . Its deviation from the thermodynamic limit is a universal function of LT. We may 
derive from Z any other thermodynamical function. The quantities we will concentrate on 
are the pressure p, the internal energy density e and the free energy density. These are given 
by, neglecting exponential corrections, 

V \ dT ) v v ; ^ 15 IV)' K J 

-/ = -logZ = (iV 2 -l)^--logv / ZTj . (5) 



According to ([T]) we have e — 3p = 0, while 



p+f _ N 2 - 1 log(LT) 3 - 1 
T 4 ~ 6 (LT) 3 



(6) 
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FIG. 1: Deviations from the thermodynamic limit (N s — > oo) of three different ideal lattice gases 
with finite spatial extent for two different values of the ratio N s /Nf. Some symbols are shifted to 
the right for clarity. For comparison, the energy density of one of these models is drawn in units 
of the continuum limit to show the dependence on the cut-off. The data are taken from Table I 
and II of Ref . [ll|] . One observes a remarkable cancellation of the UV cut-off effects for Nt > 4 and 
the emergence of a universal function of jf-. The horizontal lines are not fits, but denote the two 
values assumed by the function 1 — ^(N t /N s ) 3 as predicted by Eq. (JH). 
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FIG. 2: Deviation from the thermodynamic limit Jsb of the free energy density / of the free gluon 
gas. 



therefore we cannot trade p for — / unless LT is large enough. The presence of a log term 
in / makes it vary very slowly with the shape LT (see Fig. 2). The variation is almost 10% 
for N s /N t ~ 3 where the simulations with staggered quarks are currently performed 9| and 
reduces to 1% only at N s /N t ~ 8 which seems presently unattainable. 

A non-trivial numerical check of the above formulae comes from the mentioned study U] 
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on T — > oo limit of £77(3) gauge theory, where three different 0(a 2 ) and 0(a 4 ) improved 
actions were used to explore finite-size effects on ideal lattice gluon gases. The energy 
density e(N t /N s , N t ) was evaluated by numerical integration for some values of N t and 
for two different ratios N s /N t = 4 and N s /N t = 6 (5 significant digits) as well as in the 
thermodynamic limit e(0,N t ) (7 significant digits). Their ratios are plotted in Fig.l as 
functions of Nt- In principle, these should be model-dependent functions of N s /N t and 
N t , but there is a remarkable cancellation of the UV cut-off dependence for N t > 4. As 
a consequence, there is no other intrinsic length scale in the game besides the size of the 
system, hence these ratios are expected to collapse toward a universal function of N s /N t 
which, according to (jlj), should be 1 — ^(N/ N s ) 3 . The numerical data fit perfectly this 
prediction within numerical accuracy. 

Let us come to a proof of the main formula (J2J) . For sake of generality, I shall treat the 
case of a massless scalar field in a (hyper) cubic box of volume V = L D at a temperature 
T. The canonical partition function Z is defined by a functional integral over all periodic 
field configurations with period L in the space directions and with period (3 = 1/T in the 
imaginary time. The dimensions of the box are large compared to (3, i.e. LT ^> 1. 

There is a rich literature on a strictly related subject, the Casimir effect, where different 



methods have been developed to study various kinds of finite-size effects (see e.g. [13| for 
useful formulae and further references). It is however more instructive to account for the 
functional form of Z by exploiting some simple symmetry principles. First, we require that 
Z be dimensionless, of course. Owing to the absence of intrinsic length scales, Z should be 
a function of the unique dimensionless parameter LT. This yields invariance under the scale 
transformation 

Z = Z(LT)^ L^sL,f3^s(3 (7) 

Secondly, the D + 1 periodic box defines a cell of an infinite, regular, lattice. The physics 
should not depend on the choice of the fundamental region tiling the whole D+l-dimensional 
space by discrete translations. This requires modular invariance of the system. 

Periodic boundary conditions allow for a zero mode <fi = <p of the scalar field. In a system 
with zero modes the functional integral splits into two factors 

Z(LT) = C Q V= C [DetK]-? , (8) 

where C Q denotes the zero mode contribution while the other factor is the integral over the 



Gaussian fluctuations around the zero modes, described by the kinetic operator K. The 
latter can be written as the product over the eigenvalues A& of K in the usual form 

tDettf]-i=nAfc* , (9) 

k 

the prime here indicating that we are to exclude the zero eigenvalue. Under general grounds 
(see e.g. [l^] or [lj]], page 463) one can prove that a rescaling A^ — > s~ 2 Afc of all non- vanishing 
eigenvalues yields correspondingly [Detif]" — > s^lDetK}"^ , where J\f denotes the number 
of zero modes. In the present case K = —d 2 where d 2 is the Laplacian, and 

ak,«)=e(v) + {t) < 10) 

where the rrtj's and n run from — oo to oo. As a consequence, D = [Det — d 2 ]~^ is a function 
of L and (3 and the mentioned rescaling property of T> becomes 

V(sL,s(3) = sV(L,P) , (11) 

showing that the zero mode factorization (jSJ) spoils the scale invariance of T>. In order 
to recover the scale invariance of Z, the factor C Q should have dimensions of length. The 
only modular invariant quantity with this property can be constructed with the volume of 
the D + 1-dimensional cell, namely, (L D fi) 1 ^ ^. Choosing s = 1/L we get, aside from 
irrelevant numerical factors, 

Z= 03/L) 1/(D+1) £>(l,/3/L) . (12) 



^ion [l6,Q,[i? 



To explicitly evaluate T> it is convenient to resort to the ^-function regularization [16|,ll7|,ll8| 
of the Laplacian determinant 



D 



V(l,(3/L)=U[T,^ + (LTn) 2 ) . (13) 



m,i,n \i=l 



where the prime now indicates a regularised product. One of the virtues of the ( function 
regularization is that one can deal with regularised sums or products as they were absolutely 



convergent series and products |19l . |20L |21[ |. In (1131) we consider only the factors with n^O 
because the others generate irrelevant numerical constants. It is useful to classify the factors 
according to the number k of non vanishing m^. Denoting with the product over the set 
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{mi 7^ 0, . . . rrik 7^ 0, n ^ 0} and with pk the product over the set Ik = {mi > 0, . . . > 
0,n > 0} we have Ck = p 2 k +1 with 



J2m-+(LTny 



Ik L«=l 



and 



(14) 



logZ>(l,/?/L) = -W Wlogp* . 

fc=0 V « / 

We now use the (^-regulated product ]lm>o m ^ = ^ ^° rewrite pk in the form 

E?:>, 2 + (^n) 2 " 



p k = n n 

4_l m>0 



According to the known formula 

na 



i+ 



77r 



a 2 , e™ 



m>0 



m 2 27ra 



:i- e 



we are led to the key identity 



Pk+i =P k "e LE C Q k 



(fc) 



with 
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1 _ e ^v / S-i m ? +(LTn)2 



The (— regulated quantity 



4 fc) = ^E 



i ' 



'27rmi\ 2 f 2irn\ 2 



(15) 



(16) 



(17) 



(18) 



(19) 



(20) 



is the Casimir energy of in a box of size | x (^) fe with Dirichlet boundary conditions. 
Inserting (jl8H in f|T5|) . these Casimir energies combine to form the quantity 



D-i (D-l\ ... 

^o=E2 fc 7 4 fc) - 

fe=o \ k J 



(21) 



This is the zero-point energy of the massless scalar field in the periodic box of size (3 x L D 1 , 
which is exactly known even at finite values of LT (see e.g. 131]): 



D-l TV ' 



E n 



l d ~ 1 r 



(22) 



7T 2 
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where ((x) is the Riemann zeta function. Applying this to Eq. (j!5p and defining 

log Q = £ 2 fc+1 ( D ~ 1 ) log Q k ^ e~ 2 - LT + . . . (23) 
k=o \ k J 

we may rewrite Eq. (jl2j) in the form 

tD TV d+i \ log- 
log Z = — -^L^ D + 1) - -±2- - log Q . (24) 

This is the final result. It can be rewritten in a more evocative form 

logZ = —L E a — - £ log(l - e~^) , (25) 

k 

where Z can now be viewed as the canonical partition function of <fi in the asymmetric, 
periodic, box L D ~ l (3 in equilibrium at the "temperature" 1/L; the sum is over the mo- 
menta of the normal modes of energy uj^. The rotation, or modular transformation, of the 
fundamental D + 1 cell with respect to the standard approach to blackbody radiation has 
made it possible to highlight the finite size effects of log Z. Aside from small exponential 
corrections, it differs from the thermodynamic limit by a non-negligible logarithmic term. 
The latter is essential to enforce modular invariance, which was obvious at the beginning of 
the calculation. Its origin can be traced to the zero mode subtraction in Eq.flBJ). 

When D = 1 we recover the known partition function of a scalar massless boson on a 
2-torus. Of course, for a free gluon gas one has to multiply this result by the number of 
polarisation states and by the dimensions of the adjoint representation of the gauge group. 

In conclusion, in this paper it has been pointed out that the free energy of an ideal gluon 
gas in a periodic box is affected by non-negligible IR finite-size effects. It would be very 
interesting to observe similar effects in the interacting case. Unfortunately, in order to see 
them varying the ratio T/T c by acting on the coupling constant does not suffice: one has to 
modify the ratio N s /N t = LT. Notice that the finite-size deviation of / does not influence 
only the evaluation of the pressure but also the internal energy and the entropy. In fact 
in lattice simulations one calculates two different physical quantities: the trace anomaly 
A = e — 3p and the free energy density /, and one defines the energy density as e ~ A — 3/ 
and the entropy density as s ~ (A — 4/)/T; in view of Eq.fjSJ) these are true only in the 
thermodynamic limit. It has been observed {f| that the cut-off dependence in A is much 
smaller than for the pressure alone. This could be a clue to IR finite-size effects. 
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